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Abstract. We analyze the properties of the soliton solutions of a class of mod- 
els describing one-dimensional BEC with spin F. We describe the minimal sets 
of scattering data which determine uniquely both the corresponding potential 
of the Lax operator and its scattering matrix. Next we give several reduc- 
tions of these MNLS, derive their A^-soliton solutions and analyze the soliton 
interactions. Finally we prove an important theorem proving that if the initial 
conditions satisfy the reduction then one gets a solution of the reduced MNLS. 

1. INTRODUCTION. It is well known that Bose-Einstein condensate (BEC) 
of alkali atoms in the F = 1 hyperfine state, elongated in x direction and con- 
fined in the transverse directions y, z by purely optical means are described by a 
3-component normalized spinor wave vector <f>(x, t) = (4>i, 4>o, <I>_i) T (x, t). Con- 
sidering dimcnsionlcss units and using special choices for the scattering lengths one 
can show that &(x, t) satisfies the multicomponent nonlinear Schrodingcr (MNLS) 
equation [14], see also [15, 19, 26, 2, 22]: 

id t <S> 1 + <9 2 $i + 2(|$ 1 | 2 + 2|$ | 2 )$i + 2** = 0, 

idt<l> + d^ + 2(|$-i| 2 + |$o| 2 + |$i| 2 )$o + 2$5$!$_i = 0, (1) 

idt^-x + 3 2 $_i + 2(|$_!| 2 + 2|$ | 2 )$_! + 2$t$ 2 = 0. 

Similarly spinor BEC with F = 2 is described by a 5-component normalized spinor 
wave vector &(x,t) — ($2, $1, §0, &-2) T {x, t). For specific choices of the 
scattering lengths in dimcnsionlcss coordinates the corresponding set of equations 
for &(x,t) take the form [27]: 

*3 t $ ±2 + d xx $± 2 + 2(#, **)$ ±2 - (2$ 2 $_ 2 - 2$!$_! + $ 2 )$; 2 = 0, 
i$S±i + d xx <S> ±1 + 2(*, **)$±i + (2$ 2 $-2 - 2$,$_, + = 0, (2) 

idt*o + d xx $o + 2(5, **)$ - (2$ 2 $^2 - 2$!$_! + $ 2 )$S = 0. 
Both models have natural Lie algebraic interpretation and are related to the sym- 
metric spaces BD.I ~ SO(n + 2)/SO(n) x SO(2) with n = 3 and n = 5 respectively. 
They are integrable by means of inverse scattering transform method [4, 25, 12]. 



2000 Mathematics Subject Classification. Primary: 35Q51, 37K40; Secondary: 34K17 . 
Key words and phrases. Bosc-Einstcin condensates, Multicomponent nonlinear Schrodingcr 
equations, Soliton solutions, Soliton interactions, Reductions of MNLS. 



1 



2 



VLADIMIR S. GERDJIKOV 



Using a modification of the Zakharov-Shabat 'dressing method' we describe the 
soliton solutions [14, 17] and the effects of the reductions on them. 

Sections 2 contains the basic details on the direct and inverse scattering problems 
for the Lax operator. Section 3 is devoted to the construction of their soliton 
solutions. In Section 4 we formulate the minimal sets of scattering data L which 
determine uniquely both the scattering matrix and the potential Q(x,t). Section 

5 gives a few important examples of algebraic reductions of the MNLS. In Section 

6 we analyze the soliton interactions of the MNLS. To this end we evaluate the 
limits of the generic two-soliton solution for t — > ±oo. As a result we establish that 
the effect of the interactions on the soliton parameters is analogous to the one for 
the scalar NLS equation and consists in shifts of the 'center of mass' and shift in 
the phase. In Section 7 we prove an important theorem proving that if the initial 
conditions satisfy the reduction then one gets a solution of the reduced MNLS. 

2. The method for solving MNLS for any F. 

2.1. The Lax representation. The above MNLS equations (1) and (2) are the 
first two members of a series of MNLS equations related to the BD.I-type symmetric 
spaces. They allow Lax representation as follows [4, 10, 12] 

Ltp(x, t, A) = id x ip + U(x, t, \)tp(x, t 7 A) = 0, 

Mijj(x, t, A) = id t ip + V(x, t, \)tp(x, t, A) = 0, 



(3) 



where 



U(x,t,X) = Q(x,t) - A J, 

V(x,t,X) = V {x,t) + XV l {x,t)-\ 2 J 1 

Vx(x, t) = Q(x, t), V (x, t) = iad J 1 + ~ [ad ~j l Q, Q(x, t)] 



(4) 



For those familiar with Lie algebras I remind that, as usual, Q(x,t) and J are 
elements of the corresponding Lie algebra, which in our case is g ~ so{n + 2). The 
choice of the Cartan subalgebra element J determines the co-adjoint orbit of jj; in 
our case J is dual to ei, see [13]. It introduces grading in g = f/ ) © where 
g(o) ^ so ( n ^ xhe root system of g(°) consists of all roots of so(n + 2) which are 
orthogonal to e\\ the linear subspace is spanned by the Weyl generators E a and 
E^ a for which the roots a £ are such that their scalar products (a, ei) = 1. 
Thus the potential 

Q(x,t) = ^2 {la{x,i)E a +p a (x,t)E- a ) (5) 

a£A+ 

may be viewed as local coordinate of the above mentioned symmetric space. The 
linear operator ad jX = [J, X] and ad J 1 is well defined on the image of ad j in g. 

In what follows we will use the typical representation of so(n + 2) with n = 2r— 1 
in which Q and J take the following block-matrix structure: 

[0 f \ 
Q(x,t)=\ p s q , J = diag(l,0,...0,-1). (6) 
V p T s J 

For physical applications one uses mostly potentials satisfying the typical reduction, 
i.e. p(x,t) = q*(x,t). The vector q(x,t) for integer F — r has 2r + 1 components 

q{x,t) = ($,._!,..., $ ,--.,$-r+i) T (x,t), (7) 
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and the corresponding matrices s enter in the definition of the orthogonal algebras 
so(2r — 1); namely X G so(2r + 1) if 

2r+l / 1 \ 

X + S Q X T S Q = 0, S =Y, (-l) s+1 £ s .n+i- s , So = -so . (8) 

s =i \ i o o / 

By E sp above we mean 2r + l x 2r + l matrix with matrix elements (E sp )ij = S si S p j. 
With the definition of orthogonality used in (8) the Cartan generators Hk = E^^ — 
i?2r+2-fc.2r+2-fc are represented by diagonal matrices. 

If we make use of the typical reduction Q = (or p* = q) the generic MNLS 
type equations related to BD.I. acquire the form: 

iqt + Qxx + 1{q\q)q - (q, s q)s q* = 0. (9) 

and for r — 2 (rcsp. r = 3) coincides with the MNLS eq. (1) (resp. with eq. (3)). 
The Hamiltonians for the MNLS equations (9) are given by 

iWs = I" dx [{d x q\d x q) - (q\q) 2 + \\{f , ,s ^| 2 ) . (10) 



2.2. The Direct and the Inverse scattering problem for L. We remind some 
basic features of the scattering theory for the Lax operators L, see [10, 12]. There 
we have made use of the general theory developed in [32, 33, 29, 3, 6] and the 
references therein. The Jost solutions of L are defined by: 

lim (j){x,t,X)e iXJx = 1, lim ijj{x,t, X)e lXJx = 1 (11) 

X— > — oo x— >oo 

and the scattering matrix T(X,t) = ip~ 1 (j)(x 7 t, X). The special choice of J and 
the fact that the Jost solutions and the scattering matrix take values in the group 
SO(2r + 1) we can use the following block-matrix structure of T(X, t) 



T(X,t) 



^ ni^ —b T c x 

b+ T 22 -s B 
\ 4 B+ T s 



m 1 


B- T 




-B+ 


T 22 


s b~ 




-b +T s 


rrii 



(12) 

where b ± (X, t) and B ± (X 1 1) are 2r — 1-component vectors, T 22 (A) is 2r — 1 x 2r — 1 
block matrix, and mf(X), and cf (A) are scalar functions. Below we often use X to 
denote the matrix inverse to X. 

Remark 1. The typical reduction p(x,t) = q*(x,t) mentioned above imposes on 
T(A, ) the constraint T^(X,t) = f(X,t) for real values of A <G R, i.e. 

m+(A) = mr**(A), Si "(A) = h +'*(A), 
c+(A) = ^'*(A), S 1 +(A) = 6 1 -'*(A). 

The Lax representation (1) allows one to prove that if q(x,t) satisfies the MNLS 
(9) then the scattering matrix T(X,t) satisfies the linear evolution equation [12]: 

^-A 2 [J,T(A,i)]=0, (14) 
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or in components: 

i^-±X 2 b ± (t,X)=0, i^-±X 2 B ± (t,X)=0, 
at at 

i— L =0 1 i— 2-=0. 

eft dt 

Thus the block-diagonal matrices D ± (X) can be considered as generating functionals 
of the integrals of motion. Thus the problem of solving the MNLS eq. is based on 
the effective analysis of the mapping between the potential Q(x,t) of L and the 
scattering matrix T(X,t). 

2.3. The fundamental analytic solution and the Riemann-Hilbert prob- 
lem. The most effective method for the above mentioned analysis consists in con- 
structing the fundamental analytic solution (FAS) of L-opcrators of type (3) and 
reducing the inverse scattering problem to an equivalent Riemann-Hilbert problem 
(RHP). Skipping the details (see [11]) we just outline the construction of FAS for 
L. Obviously the FAS, like any other fundamental solutions of L must be linearly 
related to the Jost solutions. For the class of potentials Q(x,t) with vanishing 
boundary conditions there exist two FAS X ( x i t, A) which allow analytic extension 
for A £ C± respectively. For real A they are related to the Jost solutions by 

X ± (x, t, A) = <j>{x, t, X)Sf (t, A) = 1>{x, t, X)Tj (t, A)Z>± (A) , (16) 

where Tj(t,X), Df(X) and Tj(t,X) are the generalized Gauss factors of T(A,t), 
see [29, 5, 7]: 



J ' 



T(X,t) = TJD~JSj, T(X,t) = TJDjS, 

Tj(X,t) = e+^f), Sf(X,t) = e^'^), (17) 

T;±(A)=diag ((mf^m* (mf)* 1 ), 



Here 



? ± (X,t) = (T±_ ll ...,T±,...,T ± r+1 ) T (X,t), 

r-1 

= y^(r^-E ei - efe+1 + r± k E ei+ek+1 ) + r^E ei , 
fe=i 

r-l 
k=l 



(18) 



and similar expressions for (jy^iE^j . Above we have made use of the fact that 

consists of the roots {e± — e^, e\,e\+ e-kYi^X- The functions mf and n x n matrix- 
valued functions m^ are analytic for A € C±. One can check, that the analogs of 
the reflection coefficients and -r 4 are expressed by: 

B~ _ B+ . b+ . b~ 

P = — . T =—, fy=—, ^ = —- 

m 1 m 1 raj ra{ 

Remark 2. The typical reduction means that for A £ R the reflection coefficients 
are constrained by (see remark 1 above): 

p + (X) = p-'*(A), r+(A)=r-'*(A), Ael. (19) 
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There are some additional relations which ensure that both T(A) and its inverse 
T(A) belong to the orthogonal group SO(2r + 1) and that T(X)f(X) = 1. 
The FAS ^(x, t, A) are related by: 

X + (x,t,A) = x-(^,i,A)G ,j(A,0, G ,j{\,t) = Sj(\,t)S$(\,t) (20) 

Below for convenience we introduce ^(x, A) = x ± (a;, A)e lAJ:r which satisfy the 
equation 



i ^ + Q(^±(x,A)-A[J,e ± (2:,A)] =0, (21) 



and the relation 



lim ^(x,t,X) = l, (22) 

A— >oo 

Then £±(2;, A) satisfy the RHP's 

£ + (M,A)=rOM,A)G, 7 (:r,t,A), 
Gj(x,t,\) = e-' lXJ ^ + ^G ,j(X)e lXJ ^ +xt \ 

with sewing function Gj(x, t, A) uniquely determined by the Gauss factors Sf(t, A) 
taken for t = 0: 

Go,./ (A) = 57(0,A)SJ(0,A). 
The analyticity properties of these FAS follow from the equivalent set of integral 
equations: 

/x 2r-l 
dye -iK*-v) q r - p m+ +hj (y,X), 
-° P =i 



(*> A) = S kj + 1 / dtf(^_ fc+1 (y)^-(y, A) 

*/ — oo 

- (-l)^Vr+*-i(lO& +llj (»,A)), 2 < fc.j < 2r; 

2r — 1 

£+ +1>j (x,A)=£ 2r . + i,;+i [ dj/e^"^ ^(-ir+ 1 e r+p (y)^ + +1J (j/ J A), 

P =i 

(24) 

and a similar set of integral equations for 



^(x,X) = 6 lj + i f X dye-^-y) £ Qr-p(y)^p+ij(y, A), 

J — OO i 

(25) 



P =i 

£fcjO> A ) = + t / dy(q*_ k+1 (y)^ d (y, A) 

o 

2r-l 



- (-lJ^Vr+k-iGO^ufoA)), 2 < fc.j < 2r; 

e 2 - r+lj (x, A) = S 2r+1J + i fdye iX ^ £ (-ir +1 q *_ r+p (y)^(y,X), 

Jco p=i 

The RHP (23) with the additional condition (22) is known as an RHP with 
canonical normalization. 

Remark 3. An immediate consequence of the analyticity of ^(x, t, A) is that 
D ± (X) are analytic functions for A € C±. This fact follows from the relation 
)im x ^ 00 £ ± (x,X)=D ± (X). 
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Zakharov and Shabat proved a theorem [32, 33] which states that if Gj(x,X,i) 
satisfies: 

dC 

i— - \[J,G(x,X,t)] = 0, 
dx (26) 

then the corresponding solutions of the RHP allow one to construct x ± (x. A) = 
^{x, X)e~ lXJx as a fundamental solution of the Lax pair eq. (1). 

We will say that A) is a regular solution to the RHP (23) if the block- 

diagonal part of it has neither zeroes nor poles in its whole region of analyticity. 

If we have solved the RHP's and know the FAS £ + (x, t, A) then the formula 

Q(x,t)= lim x(j-£ + (x,t,X)ji + (x,t,xj), (27) 
allows us to recover the corresponding potential of L. 

3. Singular solutions of RHP and soliton solutions of MNLS. Zakharov- 
Shabat's theorem ensures that if a given RHP allows regular solution, then this 
solution is unique. However the RHP may have many singular solutions. The 
construction of such singular solutions starting from a given regular one is known 
as the dressing Zakharov- Shabat method [32, 33]. Indeed, if £y(x,t, A) are regular 
solutions to the RHP, then 

^(x,t,X) = u(x,t,X)£(x,t,X) (28) 

with conveniently chosen dressing factor u{x,t, A) may again be a solution of the 
RHP [32, 33]. Obviously this factor must be analytic (with the exception of finite 
number of singular points) in the whole complex A-plane and can explicitly be 
constructed using only the solution of the regular RHP. 

In order to obtain 7V-soliton solutions one has to apply the dressing procedure 
to the trivial solution of the RHP £o(x, t, A) = 1. We choose a dressing factor with 
2iV-poles [11]: 



N 



I + \\ n i V / A k(x,t) B k (x,t) i 
u{x,t,X) = l+y r^ + - — • (29) 

fe =i V A " A k A - A fe / 
The 7V-soliton solution itself can be generated via the following formula 

N 

Q N , s {x,t) = J2[J,A k (x,t)+B k (x,t)]. (30) 
k=l 

The dressing factor u(x, A) must satisfy the equation 
du 

*£)x + < 5 JV ' S *- X ' *) u ( a; ' *) A ) _ A I^' u ( x > *i A )] = (31) 

and the normalization condition lim,\->oo u(x, A) = 1. 
The residues of u admit the following decomposition 

A k (x,t)=X k (x,t)F k r (x,t), B k (x,t) =Y k (x,t)Gl(x,t), 

where all matrices involved are supposed to be rectangular and of maximal rank s 
[30, 9]. By comparing the coefficients before the same powers of A — Xf in (31) we 
convince ourselves that the factors F k and G k can be expressed by the fundamental 
analytic solutions x^{x,t, X) = e - 4A 0+At),/ as f vj ows 

F k T (x,t) = F k r fi l X +(x,t,X+)}- 1 , Gl(x,t) = Gl lxo(x,t,X k -)}- 1 . 
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The constant rectangular matrices F k0 and G^o obey the algebraic relations 
FkfiSoFk,o = 0, GfrfiSoGkfi = 0. 

The other two types of factors X k (x, t) and Y k (x, t) are solutions to the algebraic 
system 

l^k Ik ilk 

c n \ - XiFfSpGk v „ \ - YiGfSpGk 
i A i _A fc z^fc A * ~ A fc 

The square sxs matrices a k (x, t) and (3 k (x, t) introduced above depend on Xo an( i 
Xq and their derivatives by A as follows 

a k {x,t) = -F T fc [xo (x,t, A+)] _1 9 A X([(x,t, A^S'oFq^ +a ,fc, 

T — — 1 — — \ / 

/3fc(x,i) = -G Q j,[x (x,t,\ k )]" 9ax ( x j*) A fc )SbGo,fc + Au- 

Below for simplicity we will choose F k and to be 2r + 1-component vectors. 
Then one can show that a k — fik — which simplifies the system (32). Wc also 
introduce the following more convenient paramctrization for F k and Gk, namely 
(see eq. (35)): 



Fk{x,t) = S \rik(x,t)) 
G k (x,t) = \nt(x,t)) = 




(34) 



where i?ofc are constant 2r — 1-component polarization vectors and 

ij - v^)t -t- VOj, 



Zj = i>j(x + 2fj,jt) + £ Q j, (f>j = fijX + (/i- - vj)t + 5q 



(nJ(x,t)\So\rij(x,t))=0, or (vo,jS v OJ ) 



(35) 



With this notations the polarization vectors automatically satisfy the condition 
(n,j(x,t)\So\nj(x, t)) = 0. Thus for iV = 1 we get the system: 

\Y) = -^- X ^\ |Xl)= (A^-AD5oK) i (36) 
Wi\ni) («IK) 

which is easily solved. As a result for the one-soliton solution we get: 

qis = - 11 t (e- 2l So |^oi) +e 2l |^oi» , Ai = cosh(2*i) + <4Ki). 

(37) 
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For n — 3 we put v$k = l^ofc l eQofc and get: 

g. _ >/ 2 l^01;l t/ 01;3|(A^ - A£) j^ 1±i p la 

^ls;±l a 

^1 

x (cosh(2i =F Coi) cos(ai 3 ) - isinh(zi =p Col) sin(ai 3 )) , 

V2\v 0l , 2 \{\t - X^) _u .... (38) 

$i s:0 = — 7 e V1 (smhzi cos(a 2j + zcoshzi sm(Q 02 )) , 

&i 

Pl3 = 7,(^03 ~ a l), Coi = j. In |^"^| i "13 = ^("03 + «0l), 

Note that the 'center of mass' of $i s; i (resp. of <f>i s; _i) is shifted with respect to 
the one of $i S ;o by Coi to the right (resp to the left); besides |$i s; i| = |$i S ;-i|j i- e - 
they have the same amplitudes. 

For ?! = 5 we put v Q k = |fofc|e Qofc and get analogously: 



(39) 



$ ^ _ a/2Ki ; ii/oi ;5 |(A+-A 1 ) i<j>1±iPls 

^ls;±2 a 

^1 

x (cosh(zi =f Coi) cos(ai 5 ) -ismh(zi T Coi) sin(ai 5 )) , 

g. _ \/2|t/oi;2^oi;4|(Aj" - A£) i<t , 1±i/32i 
ls;±1 ~ Ai 

x (cosh(zi =F C02) cos(a 2 4) - isinh(zi =p Coi) sin(a 2 4)) , 

^is;0 = — V^l^oi^K^i ^ll e -i0i ( cos \ 1Zl cos(a 3) — isinhz! sin(a 3)) , 

Pis = ^(«05 - a i), Coi = ^ hi |^ 01,5 | , a>i 5 = ^(a 05 + a i), 

#24 = ^("04 - «02), C02 = j: In l^ 01 ' 4 ] , a 24 = ^(«04 + «02)- 

l I |^01;2| ^ 

Similarly the 'center of mass' of $i s; 2 and $i s; i (resp. of < &i s; _2 and $i s: _i) are 
shifted with respect to the one of $i s; o by Coi and C02 to the right (resp to the left); 
besides |$i B . 2 | = |$i s; -2| and |$ ls;1 | = |$i s; -i|. 
For N = 2 we get: 

\m(x,t)) 

\n 2 (x,t)) 
So\n*(x,t)) 
S \ri2(x,t)) 



X 2 (x,t)f 21 


ii(o:,i)«ii 


Y 2 (x,t)n 21 


x 2 - xi 


A r -A+ 


A 2 - A+ 


X 1 (x,t)f 12 


Yi(x,t)n 12 | 


Y 2 (x,t)n 22 


A+ - Aj 


X 1 - Aj 


^2 ~ ^2 


Xi (x,t)nii 


| X 2 (x,t)Ku 


, Y 2 (x,t)f^ 


X 2 - Af 


X 2 - Aj 


A 2 - A 1 


Xi (x,t)n 2 i 


! X 2 (x,t)n 22 




X 1 - A 2 


^2 ~ ^2 


X 1 - A 2 



(40) 



where 



f kj (x, t) = e **-^-*(**-*i) + e *i-** +i (**-*i) - 2 (i^so^oj) , 



(41) 
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In other words: 
/ 

fl2 



MX = 



X, —An 



A 2 — Aj 



Kl2 



V 



-AT 



^21 



A--A+ 

«12 

x ± — x 2 

o 

_/r 2 _ 

A-, — A Q 



A A2 Ag / 

Wc can rewrite .M in block-matrix form 
M21 M22 



A 2 — A + 
A 2 -Aj 













Im) \ 










|n 2 ) 










SoH) 










S \n* 2 ) J 



M = 
M n = 



/12 

A2" - 



1 

-1 




The inverse of M is given by: 

A/f 1 



M~ L = 



-N{ x MwM\ 



-Afz 1 M 2 iMi 



From eqs. (42) and (44) we obtain [12]: 
1 

Z 
1 



AC 1 



A/" 2 = M* n - M21M11M 



112 



\*i) 

\*2) 

\Yi) 
\Y 2 ) 



J12 



Z 
1 

z 
1 

z 



A x - A 2 

fl2 



K 22 c 1 *, 
A 2 - A 2 



Kl2 



A 



1 

K22 



A« 



ni + 



^21 



Aj - A 2 



2 

-|ni) 



u 

K21 



A 2 

rr|n 2 ) 



A 2 - A x 



AT 



S |n 2 



/1 



K12 



A 2 - A-l 



A x - A 2 
A+^ 



Af — Aj 



A 



|na> 



12 



where 



Z 



l/i 2 | 2 



K12K21 



A 2 - A+ 



K11K22 



(42) 



(43) 



(44) 



(45) 



ja+-a+p |A+-Arp 1 4^2 y (46) 

Inserting this result into eq. (30) we obtain the following expression for the 
2-soliton solution of the MNLS: 

Q 2s (x, t) = [J, A 1 +B 1 +A 2 + B 2 ] = \ [J, C(x, t) - S C T (x, t)S ], 



Z' 



C(x,t) = -^—\ni)(n{\ - -^I^K)<4| 

A 2 — A 2 A 2 — A j 



«21 



A 



1 _K)(4|- ,- A *\- K)(n 2 |5 - 



-A 2 " 

,/l2 



A 



AT 



AT 



at! 



a: 



\n2)(n{\ 



-SoK)KI- 



(47) 



4. The minimal sets of scattering data. It is well known that the locations of 
the sing ularities of the RHP A^ e C± are zeroes of the functions mf(X) and discrete 
eigenvalues of the Lax operator L. We will say that these eigenvalues are simple 
if the corresponding cigensubspaces are one dimensional. This corresponds to our 
choice of Fk(x,t) and Gk{x,t) as vectors. Eigensubspaces of higher multiplicities 
s > 1 can be obtained choosing Ff.(x,t) and Gk{x,t) as s x (2r + 1) matrices of 
rank s. 
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Theorem 4.1. Let the potential Q(x, t) be such that the corresponding Lax operator 
L has finite number of simple discrete eigenvalues located at the points X k G C± 
respectively, k = l,...,N. Then as minimal sets of scattering data uniquely de- 
termining both the scattering matrix T(X,t) and the corresponding potential Q(x,t) 
one can consider the sets 



Ti = {f+(A,0), AgM; rj[, A+GC+, k = l,...N}, 
% 2 = {p ± (X,0), Ael; 4, A^gC+, k = l,...N}, 
where the constant vectors fnt and " 



(48) 



Ofc 



ok ana Pok 

V2v ok , p+ k = V2p k I , (49) 



e -iok + iSok J \ g — nok+i6ok 



(50) 



and the vectors v$ k and fx^ k satisfy the normalization condition (i? Z fc s o^ofc) = 1 and 
(^ofc s oMofe) = 1- 

Remark 4. The data A^~ and A^ = (A^~ )* characterize the discrete eigenvalues of 
L. The vectors f Q k and f Q k = (fQ k )* (rcsp. p^ k and p^ k = {Po k )*) determine the 
corresponding eigenfunction of L. Note also that by definition these vectors satisfy 
(TjjfsoT^) = and (Pofe T SoPofe) = °- 

Outline of the proof. Let usbegiven71. Using r + (A, t) andr~(A,i) = (f + (X,t))* 
we construct 5^(A, t) and SQj(X,t) and therefore obtain also the sewing function 
Go(A,i) = SQ j(X,t)SQ j(X,t) for a regular RHP. According to the Zakharov-Shabat 
theorem it has unique solution £^(x,t, A). The corresponding regular potential is 
obtained by: 

Q (x, t) = lim A (j - $(x, t, X)Jit(x, t, A)) 

A— too \ / 

= [J,&(x,t)], 

where ^o\{x,t) = lim^^ \($ (x, t, A) - 1). 

Next we use the dressing method to dress the regular solution £^(x,i, A) with 
the dressing factor u(x, t, A) of the form (29). In order to do it we make use of the 
set of eigenvalues X k and X k = {X k )* and instead of the polarization vectors (34) 
we use: 

K(M)> =U(x,t,Xt)e- iX ^ x+x t t)J r ok +, 

(51) 

\K(x,t)) = ^(x,t,X k )e~^^V J T 0k -. 

After solving the algebraic equations for \X k (x 1 t)) and \Y k (x,t)} we find explicitly 
the dressed potential 

N 

Q(x,t) = Q (x,t) + Y^[J,Mx,t) + B k (x,t)], (52) 
fe=i 

which proves the first part of the theorem. 

Let us now show how one can recover T(A, t) from 71. Given the regular solution 
^ (x, t, A) we can find 

D±j{X)= lim £o(x,t, A), (53) 

' x— >oo 

and also 

(A)£>± (A) = lim e^ Xx+x *^ J £{x,t,X)e- l(Xx+x2t)J . (54) 



BEC AND SPECTRAL PROPERTIES OF MNLS 



11 



Thus we have recovered all Gauss factors T^j(X), Dqj(X) and S^j(X) of the 'un- 
dressed' scattering matrix T (X,t), so 

T ,j(X,t) = T^j(X,t)D^j(X)S^j(X,t). (55) 

In order to take into account the effect of dressing we make use of the relations 
between the dressed and undressed Jost solutions: 

ip(x,t,X) = u(x,t,X)ip (x,t, X)u+(X), 
4>{x, t,X) = u(x, t, X)cj)o(x, t, A)w_ (A), 
where u±(X) = lim. r ._ ) .-|- 00 u(x, i, A). As a result we get: 
T(A, t) = $(x,t,X)(f)(x,t,X) 

= u + (X)Mx, t, A)0o(*. t, A)w-(A) (57) 
= u+{X)T (X,t)u-{X). 
Skipping the details we state the result: 

/ c(A) \ / l/c(A) \ 

u+(A) =01 , «_(A) = 10, (58) 
\ l/c(A) / \ c(A) / 

where c(A) = l\f=i 

The fact that the set T2 is also a minimal set of scattering data is proved analo- 
gously. 

□ 

5. Reductions of MNLS. Along with the typical reduction Q = mentioned 
above one can impose additional reductions using the reduction group proposed by 
Mikhailov [21]. They are automatically compatible with the Lax representation of 
the corresponding MNLS eq. Below we make use of two types of Z2-reductions[8]: 

1) dU\x,t,X*)Cy l =U(x,t,X), CiV*(x,t,y)Ci 1 = V(x,t,\), 

2) C2U T (x, t, X)C^ 1 = —V (x, t, A), C 2 V T (x, t, X)C^ 1 = —V{x, t, A), 

(59) 

where C\ and C2 are involutions of the Lie algebra so(2r+ 1, i.e. Cf = 1. They can 
be chosen to be either diagonal (i.e., elements of the Cartan subgroup of SO(2r + l)) 
or elements of the Weyl group. 

The typical reductions of the MNLS eqs. is a class 1) reduction obtained by 
specifying C\ to be the identity automorphism of q; below we list several choices 
for C\ leading to inequivalent reductions: 

la) d = 1, p(x) = q*{x), lb) Ci = Ki, p{x) = K 01 q*(x), 

lc) d = S e2 , p{x) = K 02 q*(x), Id) C 1 = S e2 S ea , p{x) = K 03 q*(x). 

(60) 

We also make use of type 2) reductions: 

2c) C 2 = Ki, q(x) = -K 0i s q(x), p(x) = -K w s p(x), 

2f) C 2 = K 5 , q(x) = K 05 q(x), p(x) = K 05 p(x), 

where 

K,j = block-diag (1, K j, 1), K 01 = diag (ei, . . . , e r -i, 1, £r-i, ■ • ■ , d), (62) 
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for j — 1, 2, 3, 5 and 
may take the form: 



±1. The matrices A 2, K 03 and A 4 are not diagonal and 



































A02 = 


I 





-1 


i) 






K 4 = 


I 





A 04 




) 








1 















1 





, 








( 














-1 \ 




( 














-1 \ 









1 
























1 





A 02 = 










-1 








, Km = 










-1 





















1 












1 













V 


-1 











J 




V 


-1 











/ 



(63) 



Each of the above reductions impose constraints on the FAS, on the scattering 
matrix T(A) and on its Gauss factors Sf (A), Tf(X) and D^(X). For the type 1 
reductions (cases la) - Id)) these have the form: 



(5+(A*))t = A7 1 5-(A)/i, 
(£>+(A*))t=if7!D-(A)^ 



(T + (X*)y = Kj 1 T~(X)Kj 



(64) 



= K 0j T p + = K 0j p-'*, 3 = 1, 2, 3 

where the matrices Kj are specific for each choice of the automorphisms C\, see 
eqs. (60). In particular, from the last line of (64) and (61) we get: 

(m+(A*))*=m7(A), (65) 

and consequently, if mf(X) has zeroes at the points AjJ", then m 1 "(A) has zeroes at: 



k I > 



k = 1, 



and 



For the type 2) reductions we obtain: 
2c) (S ± (X)) T =K^ 1 S ± (X)K 4 
(D^f = K 4 1 D ± (X)K 4 

T ± = -A 4S f ± , 

2f) (S+(\)f = Ks 1 S-(-\)K 5 
(D+(X)f = K^D-(-X)K B 
f + (A) = A 05 r-(-A), 



(T ± (X)) T = A 4 - 1 f ± (A)A 4 
P ± = -K 0i s< i p ± , 



(T+(A)f = A7 1 T-(-A)A 5 



p + (X) = -K Q5 p-(-X), 



(66) 



(67) 



(68) 



For the 2e) reduction with n = 3 we may choose A4 to corresponds to the Weyl 
group element S ei , so A 4 = 1. As a result we get: 

$1 = -$_i (69) 

and $0 arbitrary. This reduction of eq. (1) is also important for the BEC [22]. 
From (67) we find voi = ^03- The effect of this constraint is that for the onc-soliton 
solution we get <&i s; i = — 5>is;-i- 

Our next remark following [23] is that this reduction applied to the F = 1 MNLS 
(1) leads to a 2-component MNLS which after the change of variables 



-(u?i +iw 2 ), 



V2 



(wi - iw 2 ), 



(70) 



leads to two disjoint NLS equations for w\ and w 2 respectively. 
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It is only logical that applying the constraint uqi = 1^3 the explicit expression for 
the one-soliton solution (38) simplifies and reduces to the standard soliton solutions 
of the scalar NLS. 

For the other two examples of type 2) reductions we choose n = 5 and JCj, and 
K§ correspond to the Wcyl group elements S e2 S e3 and S e2 - ea respectively. Then 
Kq4 = — So and 



( 


1 











\ 


1 






















1 






















-1 




I 








-1 





) 



(71) 



For these choices of K4, K5 we obtain: 
2e) $ 2 = $-2, 



2f) 



±2 



±- 



vr 



= 4> 



±1; 



$1 = 
$±1 = 



vT 



(72) 



It reduces the F = 2 spin BEC model into the F = 1 model. 

The corresponding relations for the Gauss factors and for the polarization vectors 
are given by: 



±2 



±- 



vr 



<i> 



:bl 



± 1 > 



(73) 



6. Two Soliton interactions. In this section we generalize the classical results of 
Zakharov and Shabat about soliton interactions [31] to the class of MNLS equations 
related to BD.I symmetric spaces. For detailed exposition see the monographs 
[29, 3]. These results were generalized for the vector nonlinear Schrodinger equation 
by Manakov [20], see also [1, 16, 28]. The Zakharov Shabat approach consisted 
in calculating the asymptotics of generic TV-soliton solution of NLS for t — > ±00 
and establishing the pure clastic character of the generic soliton interactions. By 
generic here we mean iV-soliton solution whose parameters A^ = fik ± u>k are such 
that fj,k 7^ fi>j for k 7^ j. The pure elastic character of the soliton interactions is 
demonstrated by the fact that for t —> ±00 the generic iV-soliton solution splits 
into sum of TV one soliton solutions each preserving its amplitude 2v^ and velocity 
/ifc . The only effect of the interaction consists in shifting the center of mass and the 
initial phase of the solitons. These shifts can be expressed in terms of only; for 
detailed exposition see [3]. 

Let us apply these ideas to the MNLS equations studied above. Namely we use 
the 2-soliton solution (47) derived above and calculate its asymptotics along the 
trajectory of the first soliton. To this end wc keep z\{x, t) fixed and let r = z 2 ~ z i 
tend to ±00. Therefore it will be enough to insert the asymptotic values of the 
matrix elements of M. for r — > ±00 and keep only the leading terms. For r — > 00 
that gives: 

k 2 2 — e 2r cxp(i/ 2 zi/vi) + 2Ci, 

K12 = e T exp((l + f 2 M>i + i(<l>i - <h)) + 

(74) 

K21 = e T exp((l + V2/v 1 )z 1 - i{(j>\ - 2 )) + 0(1), 
/ia = e T cxp(-(l - v % jv x )z x + i{<f>i ~ 0a)) + O(l), 
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while for t —> — oo we get: 

K22 — e~ 2T exp(— i^zi/fi) + 2Ci, 
K12 = e _T exp(-(l + vi/vx)zx - i{<t>x - <fe)) + O(l), 
Kai = e~ r exp(-(l + ^M)zi + t(^i - <fe)) + 0(1), 
/12 = e- T exp((l - y 2 /vx)zx - - <h)) + 0(l), 

After somewhat lengthy calculations we get: 

.. . , + , i v / 2^ie-^ 1 - Q +)(e^^+ So |i; 01 )+ e ^+'-+| I 7* 1 )) 

hm q2 S {x,t) = 1 . 

T ^°° cosh(2(zx + r + )) + (vl 1 ,v x) 

.. _ , + , i A /2^ie- i ^ 1 + a +)(e- J!l +'"+so|i7oi)+e 2l - r +K 1 » 
hm g 2s (a;,t) = 7 , 

cosh(2(z! - r+)) + 

where 

A+ - A 



(75) 



(76) 



hi 



A l ~~ A 2 



A : - A 2 



arg 



A 1 - A 2 

For n — 3 and n = 5 the right hand sides of (76) coincide with the one-soliton 
solutions (38) and (39) respectively. This means that the 2-soliton interaction for 
the above MNLS eqs. is purely elastic. The solitons preserve their shapes and 
velocities and the only effect of the interaction consist in shifts of the center of mass 
and the phase. From this point of view the interaction is the same like for the scalar 
NLS eq. 

It is important to check whether the A-soliton interactions consist of sequence 
of elementary 2-soliton interactions and the shifts are additive. 

7. Effects of reductions and initial conditions on MNLS. 

Theorem 7.1. Let the minimal set of scattering data Tj, 3 = 1,2 for t = satisfy 
the reduction conditions (67). Then the solution q{x,t) of the MNLS with such 
initial data will satisfy the corresponding reduction 2e) (61). 

Proof. Let the minimal sets of scattering data, say Tx satisfy the reduction condi- 
tions (67) for t = 0. It is easy to check that their evolution law (15) is compatible 
with the reduction, so (67) will hold for all t > 0. As a result the corresponding 
Gauss factors S , T and D^, and consequently, the sewing function in the RHP 
G(x,t,X) will satisfy 

G{x,t, A) = K^ l G T {x,t,X)K i . (77) 

The next consequence is that both ^ and K^ 1 ^* 1 - A4 are solutions of the RHP (23) 
with the same sewing function and the same canonical normalization. Therefore 
from the uniqueness of the solution of RHP we get that the regular solutions of this 
RHP satisfy: 

&(x, t, A) = K^' T (x, t, A)A 4 . (78) 

Next we note that the scattering data related to the discrete spectrum also satisfy 
the reduction conditions. This means that the dressing factor u(x, t, A) and the 
singular solutions £^(x, t, A) = u(x,t,X)^(x,t,\)u-(\) also satisfy: 

u(x,t,X) = A7 1 u T (a;,t,A)A:4, ^(x, t, A) = K^ l ^ ±%T (x, t, A)A 4 . (79) 

It remains to check that from equations (27) and (79) there follows: 

Q(x, t) = -K^Q T (x, t)K 4 . (80) 
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Remark 5. Note that the above arguments are not specific for the choice of K4. 
The above theorem can be proved along the same lines for any reduction of type 1 
and type 2. 

A simple consequence of the above theorem is the following. Consider n = 3 and 
choose Tj + = 73 for t = 0. Then the corresponding solution of F = 1 BEC (1) will 
also satisfy $1 = — $_i for all t > 0, i.e. will be a solution to 

idt&i + d 2 x <Z>! + 2(|$!| 2 + 2|$ | 2 )$! - 2$^$g = 0, 

idt*o + + 2(2|$x| 2 + |$ | 2 )$o - 2$*$ 2 = 0, 

If we insert eq. (70) into (81) we obtain 

idtWi + diw\ + 2\wi\ 2 wi = 0, 

(82) 

id t w 2 + dlw 2 + 2\w 2 \ 2 w 2 = 0, 

Therefore, if we want to analyze the specific features of F = 1 BEC we have to 
avoid such initial conditions. 

Similarly, if for n = 5 we choose in (39) foi-x = foi,5i ^oi:2 = — "01,4 w e will 
obtain in fact a solution to F = 1 BEC. 



8. Conclusions and discussion. Using the Zakharov-Shabat dressing method 
we have obtained the two-soliton solution and have used it to analyze the soliton 
interactions of the MNLS equation. The conclusion is that after the interactions the 
solitons recover their polarization vectors vok, velocities and frequency velocities. 
The effect of the interaction is, like in for the scalar NLS equation, shift of the 
center of mass z\ — > z\ + r + and shift of the phase <j)\ — > cf>i + a + . Both shifts are 
expressed through the related eigenvalues only. 

The next step would be to analyze multi-soliton interactions. Our hypothesis is 
that each soliton will acquire a total shift of the center of mass that is sum of all 
elementary shifts from each two soliton interactions. Similar result is expected for 
the total phase shift of the soliton. 

Finally we have proved a theorem, stating that a symmetry imposed on the 
minimal set of scattering data leads to a symmetry of the corresponding solution. 
So if we want to analyze the specific features of a given MNLS we have to avoid 
such initial conditions. 
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